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Abstract 



O . 

^ . Quasifree 77 photoproduction from nuclei is studied in the Distorted Wave 

Impulse Approximation (DWIA). The elementary eta production operator 
contains Born terms, vector meson and nucleon resonance contributions and 
^ ' provides an excellent description of the recent low energy Mainz measurements 

■ on the nucleon. The resonance sector includes the 5'ii(1535), i-*ii(1440) and 

1)13(1520) states whose couplings are fixed by independent electromagnetic 
and hadronic data. Different models for the r]N t-matrix are used to construct 
a simple rjA optical potential based on a tp-approximation. We find that the 
exclusive ^(7, rjN)B process can be used to study medium modifications of 
the A^* resonances, particularly if the photon asymmetry can be measured. 
The inclusive ^(7,r/)X reaction is compared to new data obtained on ^-^C, 
^^Ca, and is found to provide a clear distinction between different models for 
the r]N t-matrix. 
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I. INTRODUCTION 



Photoproduction of the isoscalar eta meson with mass of approximately 550 MeV offers 
the possibihty of investigating the rjN- and ?774-interaction, as well as the basic production 
amplitude. Recently it has been suggested by Wilkin that the very large production cross 
sections found in the pd — > rfHe reaction near threshold could be due to an rjN scattering 
length that is much larger than the value extracted from a coupled-channels analysis of the 
TT~p — > rjn and vrA^ — > vrA^ data 0. Based on the K- matrix formalism, Wilkin was able 
to reproduce the strong threshold energy dependence of pd TfHe by assuming an rjN 
scattering length with a real part more than twice as large as in previous analyses. Should 
this conclusion turn out to be true, it would have dramatic implications for the r^-nucleus 
interaction. Most importantly, a larger value for the scattering length might indicate a 
larger probability for the presence of a "bound" ?7-state for lighter nuclei such as ^He than 
had been expected. In fact, this seems to be required to explain the cusp-like structure seen 
for near-threshold production in the pd — ^ ^HeX reaction for a missing mass close to the 
?7-mass. This effect may also remove the discrepancy between experiment and theory for the 
reaction ^He{Tc~ ,riyH |@. 

Over the last several years there has been renewed interest in the production of r^-mesons 
with protons, pions and electrons and their interaction with nucleons and nuclei. One of 
the first r^-nuclear experiments, performed at SATURNE in 1988 [Q, reported surprisingly 
large eta production rates near threshold in the reaction p{d,^ He)ri. These large cross 
sections permitted not only a more precise determination of the 77-mass |^ but were also 
used to perform rare decay measurements of the eta 0. Additional experiments involving 
pion induced eta production were performed at Los Alamos [0]. Again, the experimental 
cross sections at threshold region of the reaction ^He(7r~, ?7)'^H are above the theoretical 
calculations 

The advent of high duty-cycle electron accelerators opens for the first time the op- 
portunity to study the reactions N{'j,ri)N and N{e,e'ri)N in greater detail. Our present 
knowledge of the (7,77) process is based solely on some old measurements around 20 years 
ago p], along with very few more recent data from Bates and Tokyo |Ty]. Over the last 
two years eta photoproduction from the nucleon has been measured at Mainz and at Bonn 
with an accuracy of more than an order of magnitude better than in older experiments. At 
Mainz, the TAPS collaboration has obtained high-quality data for angular distributions and 
total cross sections for photon energies between threshold and 790 MeV both on the nu- 
cleon [0] and on nuclear targets ^^C, ^°Ca [jl2[. That may be considered to be a qualitative 



breakthrough in the experimental field. Similarly, data at the higher energies up to 1150 
MeV will be provided soon by the Phoenics collaboration at ELSA |jT3[ . 

In this work we present a study of the quasifree eta photoproduction on nuclei ^(7, riN)B 
(where B is some discrete nuclear state) in the Distorted Wave Impulse Approximation 
(DWIA). The model allows for the study of the production process and the final-state 
interactions without being obscured by the details of nuclear transition densities, unlike the 
reaction of the type ^(7, ri)B where nuclear transition densities play an important role. This 
is mainly due to the quasifree nature of the reaction: the momentum transfer to the residual 
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nucleus can be made small. The key ingredients of the model are a) the single particle wave 
functions and spectroscopic factors, b) the elementary eta photoproduction amplitudes, and 
c) the final state interactions. Information on each ingredient can be taken directly from 
independent studies. Previous studies of such quasifree reactions on the pion have proven 
to be successful [|14|-[T6|. Other studies on eta photoproduction from nuclei can be found, 
for example, in Refs. [Il7|-p3[]. 

In Sec. m we briefly describe the [^,7]) operator and show comparisons with new data 
from Mainz. The DWIA formalism for quasifree eta photoproduction on nuclei is derived 
in Sec. |ITT[ In Sec. |V| we present predictions for differential cross sections of the exclusive 
reaction A{'j,7]N)B. In Sec. |V| we compare our calculations with new data from Mainz on 
the inclusive reaction A{'-f,ri)X. In Sec. ^ we summarize our findings and present a brief 
outlook. The transformation of CGLN amplitudes from the cm. frame to an arbitrary 
reference frame is given in the Appendix. The transformation is needed to implement the 
elementary amplitudes in nuclear calculations. 



II. ETA PHOTOPRODUCTION ON THE NUCLEON 



Nucleon resonance excitation is the dominant reaction process in (7,//). In contrast to 
pions which will excite A(T = 3/2) as well as N*{T = 1/2) resonances, the t] meson will 
only appear in the decay of N* resonances with T = 1/2. In the low-energy region this 
is dominantly the S'ii(1535) state that decays 45-55% into rjN, the only nucleon resonance 
with such a strong branching ratio in the t] channel. This result is even more surprising as a 
nearby resonance of similar structure, the S'ii(1650) has a branching ratio of only 1.5%. The 
other resonances that play roles in the low energy region are the Pii(1440) and Di3(1520). 

Most attempts to describe eta photoproduction on the nucleon l2^j2^ have involved 
Breit-Wigner functions for the resonances and either phenomenology or an effective La- 
grangian approach to model the background. These models which contain a large number of 
free parameters were then adjusted to reproduce the few available data. In a very different 



approach, Ref. derived a dynamical model which employs nN ttN, tcN — > irnN and 
TT~p rjn to fix the hadronic vertex as well as the propagators and the —>■ ttN to 
construct the electromagnetic vertex. This calculation represents a prediction rather than 
a fit to the jN rjN reaction. The model was later extended to include the background 
from s,M-channel nucleon Born terms and p,uj exchange in the t-channel ||2^. Due to the 



fact that 1] is an isoscalar meson, two coupling schemes are possible at the riNN vertex: 
pseudovector(PV) and pseudoscalar(PS). The latter one is not ruled out by LET as in the 
case of (7, vr). Since the resonance sector is fixed in this approach and also the vector meson 
couplings can be obtained from independent sources, one can use this model to determine 
the nature of rjNN coupling: both the coupling scheme and the coupling constant. 

In contrast to the vrA^-interaction, little is known about the r/iV-interaction and, con- 
sequently, about the rjNN vertex. The uncertainty regarding the structure of the rjNN 
vertex extends to the magnitude of the coupling constant. This coupling constant g^^^/An 
varies between and 7 with the large couplings arising from fits of one boson exchange 
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potentials [|28|. /^From SU(3) flavor symmetry, all coupling constants between the meson 
octet and the baryon octet are determined by one free parameter a, giving 

^ = 1(3_4«)2^, (1) 
An 3 An 

resulting in values for the coupling constant between 0.8 and 1.9 for commonly used values 
of a between 0.6 — 0.65, depending on the F and D strengths chosen as the two types of 
SU(3) octet meson-baryon couplings. Other determinations of the rjNN coupling employ 
reactions involving the eta, such as n^p rjn, and range from 0.6 - 1.7 ^9]. Smaller values 



are supported by A^A^ forward dispersion relations [^] with g^^j^/An + g^,j^j^/An < 1.0. 



There is some rather indirect evidence that also favors a small value for gn^N- In Ref. [pi 
Piekarewicz calculated the n-r] mixing amplitude in the hadronic model where the mixing 
was generated by NN loops and thus driven by the proton-neutron mass difference. To 
be in agreement with results from chiral perturbation theory the riNN coupling had to be 
constrained to the range 0.32 - 0.53. In a very different approach, Hatsuda evaluated 
the proton matrix element of the flavor singlet axial current in the large Nc chiral dynamics 
with an effective Lagrangian that included the ?7a(1) anomaly. In this framework, the EMC 
data on the polarized proton structure function (which have been used to determine the 
"strangeness content" of the proton) can be related to the rj'NN and the rjNN coupling 
constants. Again, his analysis prefers small values for both coupling constants. Nevertheless, 
from the above discussion it seems clear that the rjNN coupling constant is much smaller 
than the corresponding nNN value of around 14. 

In Fig. |l], we show the sensitivity of the total cross section close to threshold when varying 
the coupling constant from to 3 for the PS and from to 10 for the PV form. There is 
a large variation of more than a factor of 2 at 750 MeV for the PS case while changing the 
value with PV structure modifles the cross section only by a relatively small amount. 

Here the very precise results of the new Mainz experiment can clearly distinguish between 
the different models preferring a pseudoscalar coupling scheme with a best-flt couphng con- 
stant of g'^^j^/ An=O.A. This is further supported by comparing with the differential cross 
section as shown in Fig. |^. There is a clear distinction in the forward-backward asym- 
metry of the angular distribution between the PS- and PV-model, and PS coupling with 
g'^j^j^/ An=O.A is strongly supported. The variation in the angular distributions is due to 
the p-wave multipoles. In particular, the Mi_ multipole changes sign between PS and PV 
coupling. 

Fig. 1^ shows the complete view of the differential cross section and the photon asymmetry 
of 7] photoproduction on the proton from threshold up to 800 MeV for PS coupling with 
g'^j^j^/An=O.A. The resonance dominance near threshold is clearly seen in the cross section. 
The angular distribution is relatively flat in the cross section but peaked in the photon 
asymmetry. 



III. DWIA MODEL FOR r] PHOTOPRODUCTION FROM NUCLEI 
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A. Kinematics 



In this section we define our coordinate system and discuss some kinematic aspects of 
tlie exclusive reaction (7, r]N) and tlie inclusive reaction (7, t]) on nuclei. In the laboratory 
frame, the four-momenta of the incoming photon is k'^ = [E^, k) while the outgoing eta 
and nucleon have four- momenta of q'^ = {Er,,q) and = {E^,p) respectively. The target 
nucleus is at rest with mass Mi and the recoiling residual nucleus of mass Mf has momentum 

Q = k-q-p (2) 

and kinetic energy Tq = ^jj- The momentum transfer Q is sometimes also called the 
missing momentum. Overall energy conservation requires 

E., + Mi = E^ + EN + Mf + Tq. (3) 



As shown in Fig. ^, the z-axis is defined by the photon direction k, and we choose the 
azimuthal angle of the eta, 0^ = 0, by defining y = k x q and x = y x z. 

We assume that the reaction takes place on a single bound nucleon with four -momentum 
p'^i = {Ei, Pi) and that energy and momentum are conserved at this vertex (i.e., the impulse 
approximation). Thus p, = — Q and Ei = E^j + E^q — E^. If Q does not vanish, the struck 
nucleon is off its mass shell. For exclusive reactions this is the only reasonable off-shell choice 
since the photon, the eta, the outgoing nucleon, and the recoiling nucleus are all external 
lines and must be on their respective mass shells. For the inclusive reaction we make the 
same choice. 

The magnitude of the momentum transfer to the recoiling nucleus has a wide range, 
including zero, depending on the directions of the outgoing eta and nucleon with respect 
to the incident photon beam. However, since the reaction amplitude is proportional to the 
Fourier transform of the bound state single particle wavefunction, it becomes quite small for 
momentum transfers greater than about 300 MeV/c. Thus for all but the lightest nuclei we 
can safely neglect the nuclear recoil velocity and generate optical potentials for the outgoing 
particles in the laboratory frame. As noted in |T^, reactions of the form {'y,riN) on nuclei 



offer great kinematic flexibility and by appropriate choices one can investigate the production 
operator, the bound state wavefunction, or the flnal state interaction of the outgoing meson 
and nucleon. 



B. The DWIA cross sections 



In this section we give the DWIA formalism for calculating the quasifree eta photopro- 
duction from nuclei. Following standard procedure, the differential cross section is given 

by 
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where X) means sums over final spins and average over initial spins. The 5-function ensures 
overall energy-momentum conservation. In impulse approximation, the many-body matrix 
element can be written as a sum over single particle states a = {nljm} 

a 



where is a destruction operator and i-y^ the eta photoproduction on-body operator. The 

overlap Sa = {Jf,mf\aa\Ji,miy is conventionally called the spectroscopic factor which can 
be determined from electron scattering. After integrating over dQ and performing the sums, 
one arrives at the differential cross section 

where the single particle matrix element T = {r]; N\t^^\a;j) is given by 

T{a,X,m,) = I dV*W(r,-p) ^^^r, -q) t^,(A, k, p,, q, p) *„(r) e^"^-". (7) 



In this equation A is the photon polarization, the spin projection of the outgoing nucleon, 
and (f>!^^ the distorted wavefunctions with outgoing boundary condition, and ^q, the 
bound nucleon wavefunction. 

Another useful observable, the photon asymmetry, is defined by 



where _L and || denote the perpendicular and parallel photon polarizations relative to the pro- 
duction plane (x-z plane). We have used the short-hand notation d^a = d^a / dEj^ dflndflN ■ 



C. The nucleon single particle wavefunction 



As discussed in the section on kinematics, we are primarily interested in cases of low 
momentum transfer to the recoiling nucleus. Thus our choice of single particle wavefunctions 
for the bound state is not critical as long as the basic size of the orbital is described correctly. 
For convenience we will use harmonic oscillator wavefunctions which have the advantage 
that their Fourier transforms are simply obtained. For each nucleus under consideration, 
we adjust the harmonic oscillator range parameter until the rms radius of the ground state 
charge distribution agrees with the experimentally determined values. 

For the continuum nucleon wavefunctions we solve the Schrodinger equation with an 
optical potential present whenever the particle is to be detected. For exclusive reactins, we 
use the experimentally determined separation energies for a given orbital in order to fix the 
value of the mass of the recoiling nucleus Mf. When we are studying the inclusive (7,77) 
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reaction we allow for all possible final states by using a plane wave for the proton and taking 
the mass of the recoiling nucleus to be equal to the ground state of the A-1 system. In 
previous work on(e, e') reactions from nuclei in the quasielastic region |^3[ we found that an 
assumption similar to this described the experimental results very well. 

Many optical models for the outgoing nucleon are available. We use a non-relativistic 
reduction of the global optical model in Ref. ||3^. This model has the advantage that it fits 
nucleon scattering over a wide range of energy and A values, and hence is very useful for 
making surveys of a wide range of possible experimental situations. Once experimental data 
is available for the exclusive reaction, an optical model specific to the nucleus and energy 
range of the outgoing nucleon can be substituted. 



D. The r] optical potential 



In order to reduce the r;- nucleus many body problem to an equivalent potential scattering 
problem, we construct a simple optical potential. Distorted t] wavefunctions can be obtained 
by solving the Klein-Gordon equation 

[_V2 + _ ^2]^^^^) ^ -2cu,\/(r)0,(r) (9) 

where fi is the rj reduced mass and tu^ its total energy. We choose a simple tp approximation 
to construct the potential 

- 2urjV{r) = bp{r) (10) 

where p is the nuclear local density. This approach is justified in the low energy regime since 
S- waves dominate via the S'ii(1535) state and P-wave and D-wave contributions are very 
small. The parameter b is related to the t]N t]N scattering amplitude by 

b = 4T,Pl^f (11) 

Pcm 

where p denotes the rjN two body momentum in the respective frame. Here we consider two 
models for the r/N scattering amplitude. 



The first model is from the coupled channel approach in Ref. ||2^. We can extract the 
rjN t-matrix (in the riN cm. frame) by 



/ = ^J2/3 [ts,, + U,, cose + 2tn,,P2{^ose)] /ql (12) 

where 6 is the cm. angle and 

= [{W^ - (m + m^y){W^ - (m - m^)^]^/y2W. (13) 

The partial wave scattering amplitudes are given by 
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i.^^W) = -^-T^^{W)D{W), (14) 



where D(W) is the resonance propagator (see Eq. 2 in Ref. ||2^) and Fj^ is related to the rj 



self-energy (see Eq. 3 in Ref. |^) by = — 2ImS^. We call the optical potential based on 



this amplitude DWl. 

The second model is given in Ref. ||3^. In this analysis the nN r]N and t]N r]N t- 



matrices are obtained in a unitary, coupled, three-channel approach with the third channel 
ttttN being an effective two-body channel which represents all remaining processes. The 
7iN elastic phase shifts and the weighted data base of the ttN — >■ r/N total and differential 
cross sections are chosen as the input for the fitting procedure. The resulting rjN t-matrix 
describes the data fairly well. For numerical calculations, one can use the following scattering 
length and effective range expansionwhich accurately represents the original dominant 6*11 
amplitude: 

/ = — — - (15) 

with 

= (0.876 + ^0.274) fm, = (-1.682 - i 0.139) fm. (16) 

We call the optical potential based on this amplitude DW2. It has been pointed out in several 
references [||^ that a fit to total 7r~p — > rjn data based on the optical theorem requires 
the imaginary part of the r]N scattering length to be larger than roughly 0.22 fm. This 
would be difficult to reconcile with DWl which has /ma^=0.16 fm or the model of Ref. 
with Jmar,=0.19 fm. The latter two models were obtained by using vr^p — ^ rjn differential 



cross section data, some of those have been critized for experimental problems [35]. A new 



TT p rjn experiment at Brookhaven [BBl has recently remeasured both differential and 



total cross sections at threshold and should clear up this controversy. 

Fig. 1^ shows the energy dependence of the real and imaginary parts of the parameter b, 
defined in terms of / in Eq.(ll), for the two models. We see that the real parts are roughly the 
same for energies larger than 50 MeV, but the imaginary parts, which determine absorption, 
are quite different. As we will see later, it is the imaginary part that most strongly influences 
inclusive reactions from nuclei. It is interesting to note that for small eta energies the real 
part of b is positive (more so for DW2 than for DWl) which indicates an attractive potential 
(see Eq. (p!OD). This has led to the suggestion that the 17- nucleus system may form hadronic 
bound states [^7H 39| . In the following investigations, unless otherwise noted, we will use the 



DWl ?7-nucleus optical model. 



IV. THE EXCLUSIVE REACTION A{j,r]N)B 

The exclusive reaction A{'~f,riN)B where B is in some specific final state offers a wide 
range of experimental possibilities. In order not to be overwhelmed with details, we will 
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concentrate on cases where the momentum transfer to the recoihng nucleus is allowed to 
vary freely. This gives the widest kinematic phase space and facilitates discussions on the 
energy and angular distributions. Within our model, the reaction occurs on a single nucleon 
which is in some specified orbital in the nucleus. The nucleon can be either a proton or 
a neutron and it has some momentum distribution given by the Fourier transform of its 
wavefunction. For simplicity, we consider the coplanar setup {i.e., (pn = 0, (j)N = 180°) 
which in general leads to larger cross sections than out-of-plane setups. 

In Fig. 1^, the effects of final-state //-nucleus and p-nucleus interactions from different 
nuclei are shown for the exclusive cross section and the photon asymmetry. The results 
correspond to knocking out a proton in the P3/2 orbital of the target nuclei whose momentum 
distribution is largely responsible for the double peaks in the cross sections. We see more r) 
distortions at lower eta energies than higher eta energies, while the opposite is true for proton 
distortions. For fixed photon energy, the proton energy decreases as eta energy increases, 
see Eq. (Q). As expected, there are larger distortions in "^^Ca than in ^^C. It is interesting 
to note that distortions have no effects on the photon asymmetry, which is the same result 
found earlier in investigating exclusive (7,7rA^) reactions from nuclei [0]. Furthermore, the 
photon asymmetry is almost independent of the target it is produced from. This makes it 
a sensitive observable for studying the production process in the nuclear medium without 
being obscured by distortion effects and overall normalizations. 

In Fig. 1^ the effects of the complete and just the imaginary part of ry-nucleus final-state 
interaction are shown as a function of the eta scattering angle for eta mesons with 30MeV 
kinetic energy. Even for this exclusive reaction, most of the effects on the cross section are 
due to the absorptive part of the scattering amplitude. For inclusive reactions, where the 
relative angle between the eta and nucleon is not so important, the effects of the real part 
of the optical potential will be even smaller. This insensitivity to the real part of the optical 
potential can be used to determine the imaginary parts of different ?7-nucleus potentials. All 
of these have no effects on the photon asymmetry as evidenced in Fig. |^ and as discussed 
above. 

In Fig. H we show the sensitivity to possible medium modifications in the cross section 
and photon asymmetry for ^^C(7, rjpY^Bg s. resulting from creating an eta from a ^3/2 proton 
orbital. The results are presented as a function of photon energy at fixed eta energy and fixed 
eta and proton angles. Under this situation, the proton energy increases and the momentum 
transfer decreases as photon energy increases. Since we are interested in the relative changes, 
the results were calculated using plane waves. The solid line is the standard calculation, 
while the dashed and dotted line have the mass of the 5*11(1535) increased and decreased 
by 3%, respectively. We see large effects in the cross section and some effects in the photon 
asymmetry. 

Fig. P shows the effects of varying the mass of the Z^is resonance. We see large effects in 
the photon asymmetry while almost no effects in the cross section. In particular, the photon 
asymmetry shows a large sensitivity to a decrease in the mass of the D13. For the bare 
process this sensitivity of the photon asymmetry to the D13 resonance was first reported in 
Ref. [B^. We suspect some interplay between neighboring resonances to be involved. On 



the other hand, we found no sensitivity to varying the mass of the Pn resonance. Clearly, 
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the measurement of the photon asymmetry for the basic reaction and from nuclei would be 
very valuable in looking for possible medium modifications. 



V. THE INCLUSIVE REACTION ^(7, r?)X 



Recently, new and good quality data have been obtained at Mainz for the inclusive 
reaction on nuclear targets ^^C, ^^Ca |12|. In this section we will compare our model with 



the data. 

In the inclusive reaction the final nucleon is not observed, thus leading to all possible 
final nuclear states. As mentioned earlier, we fix the kinematics for the inclusive reaction 
by assuming the final residual A-1 nucleus is in its ground state, and we use plane waves 
for the outgoing unobserved nucleon. In our quasifree model of the reaction, the inclusive 
reaction cross section is obtained by integrating out the nucleon solid angle dflN in Eq. (H) 
and summing over all nucleon states in the target 

d^(^ = V f / "^'^ ^ (17) 

dTr,dQ^ {r{\J dE^dn^dVtN ^ ) ■' ^ ' 

Note here that the eta production amplitude is different on the proton and the neutron. In 
performing the sum over neutron states, we use the ratio Onjo^ = 2/3 extracted from the 
Mainz data for eta photoproduction on the deuteron [^]. 

The kinematics in this case can be solved by specifying the photon energy E^^ the eta 
kinetic energy and the eta angle We examined the distribution of cross sections as a 
function of eta angle and energy at fixed photon energy, and found that most of the cross 
sections are concentrated in the forward directions (^^ < 30°), unlike pion photoproduction. 
It is due to the relatively large mass of the eta as compared to the pion. At larger angles, 
cross sections are quickly reduced because of the large momentum transfer involved. Cross 
sections as large as around 0.2 yhjMeV sr^ can result at eta energies around 100 MeV. Also, 
etas with energies higher than a certain value (depending on the target) are kinematically 
forbidden. 



In Figs. ITO| and O we investigate the effect of the two different final state interactions 



on two different inclusive cross sections. In Fig. ^ we show the effects of the final-state 



//-nucleus interaction as a function of eta energy at fixed photon energy and eta angle. 



while in Fig. ^ as a function of eta angle at fixed photon energy and eta energy. The two 
different final state interactions investigated are based on models DWl and DW2 discussed 
in Section III.D. The primary result in both cases is a reduction of the cross section due to 
the absorptive part of t he optical potential. As noted earlier, the optical potential based 
on DW2 is considerably more absorptive than the one based on DWl, which is refiected in 
these figures. 

The integrated cross section da/dTr^ is obtained after a further integration over the eta 
solid angle dVL^^ in Eq. (|1^. The kinematics are solved by specifying two variables: E^ 
and T^. Fig. O shows the comparison of our calculations with the Mainz data on ^^C and 
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Ca as a function of the eta energy at fixed photon energy. The plane wave calculation 
clearly over-predicts the data by about a factor of 2, while the distorted wave calculation 
using DWl agrees with the data reasonably well, and the distorted wave calculation using 
DW2 underestimates the data by about 20 to 30%. The shape of the energy distribution is 
reproduced very well for and reasonably well for ^^Ca. 

The integrated cross section da/dO^jis obtained by integrating over the eta energy dT^i 
in Eq. ([T7|) and multiplying by the angular factor 27rsin0. The kinematics are solved by 
specifying two variables: and 6'^. Fig. |1^ shows the angular dependence at fixed photon 
energy for the comparison of our calculations with the Mainz data on ^^C and ^^Ca. While 
the agreement with is also very good here, we predict a somewhat different shape for 
"^^Ca than is measured. At lower scattering angles, the DW2 optical model furnishes a better 
description of the data, while at larger angles, the DWl optical model is much better. It 
could be that the real parts of the optical potentials are playing some role here. 

To obtain the total cross section in our model for the quasifree eta photoproduction 
from nuclei, we need to perform extensive integrations which are 4-fold altogether (the d(j)r] 
integration gives 2-7r due to the azimuthal symmetry after dfl^ is performed). These inte- 
grations are done numerically using Gauss' method. The cross sections are concentrated 
in a limited region of phase space, outside of which they are suppressed due to large mo- 
mentum transfers. We adjusted the integration limits, the number of integration points, 
and the number of partial waves, until convergence in the total cross section is achieved. 
We arranged the program in such a way that the differential cross sections discussed above 
d'^a/dTrjdflj^, da/dT^^, da/dOrj, are saved as intermediate results in the same run for the total 
cross section. We calculated at three photon energies -E^=700 MeV, 750 MeV, 780 MeV 
and interpolated to other energies since the total cross section is expected to be a smooth 
function of E^. The results are shown in Fig. [1^ along with experimental data from Mainz. 
We see that the model does a good job in reproducing the data. The distortion (absorption) 
of etas in the final state is clearly needed to achieve this agreement. Again, the absorption 
present in DWl is in better agreement with the data than the absorption in DW2. The 
absorption present in DW2, calculated in our model, appears too strong. 



VI. CONCLUSION 



We have developed a DWIA model for quasifree eta photoproduction on nuclei 
A{'y,T]N)B. The three key ingredients of the model are single particle bound state wave 
functions; the elementary eta photoproduction amplitudes; and the final state interactions. 
All of these ingredients enter the calculation in a physically transparent way so that one 
can use different models for each ingredient from independent studies. Furthermore by us- 
ing the kinematic flexibility present in the reaction, one can emphasize one or more of the 
ingredients. 

In particular, with the exclusive reaction, kinematics can be chosen so that the mo- 
mentum transfer to the recoiling A-1 nucleus is quite small and hence the calculated cross 
sections are quite insensitive to the bound state wavefunction. We found that the exclusive 
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reaction cross sections and, in particular, the photon asymmetry are very sensitive to pos- 
sible medium modifications of the intermediate nucleon resonances. We recommend that 
exclusive cross sections be measured, and once polarized photons are avaiable in this energy 
range, that the photon asymmetry be measured. This reaction seems to be the cleanest one 
available for determining whether or not intermediate resonant states are different in the 
nuclear medium than in free space. 

While exclusive cross section data is not yet available, there has been a series of inclusive 
cross section measurements from Mainz. We have compared our model to the recent inclusive 
data on and ^°Ca from Mainz and we find good agreement in our model when we use 
an ?7-nucleus optical potential based on the work of Bennhold and Tanabe [EBl. A second 



optical potential developed in Ref. appears to have too much absorption. We conclude 



that r] photoproduction from nuclei provides a good tool to extract properties of the rjA 
interaction. 
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APPENDIX A: TRANSFORMATION OF CGLN AMPLITUDES TO AN 
ARBITRARY REFERENCE FRAME 



Consider a general process of meson photoproduction on the nucleon in which the final 
baryon mass m2 can be different from the initial mass mi. Energy- momentum conservation 
gives k+pi = q+p2- The kinematics can be described by the Lorentz-invariant Mandelstam 
variables s = {k+piY, t = [k — q)"^, u = (A;— ^2)^ with the constraint s+t+u = m^+ml+m^, 
or by the invariant mass W = \/s and the scattering angle 9. 

In the center-of-mass (CMS) frame we denote the four-momenta as = (i?^,k), pi = 
{El, — k), q = {Eg, q), P2 = {E2, — q)- These CMS variables can be expressed in terms of W: 
El = (W^ + ml)/{2W), E2 = (W^ + ml- ml)/{2W), etc. The cross section in the CMS 
frame can be written as 

da/dVL = {k/q)\Tf,\'^ (Al) 

where we have used A; = |k| and q = |q| to denote the initial and final 3- momentum in the 
system, respectively. 

The transition matrix element can be expressed as 

Tf^ = {x2\F\xi) (A2) 

where xi, X2 are Pauli spinors and 

F = icr ■ e Fi + cr ■ q a • (k X e) F2 + ia • k 0- ■ e F3 + icr ■ q q ■ e F4 (A3) 



The F's are the CGLN amplitudes Note that there is a factor 17111712/ {'^'kW) absorbed 
in the F's under the normalization of Eq. ([Al|). 

On the other hand, the transition matrix element can be expressed in the invariant form 



Tf, = Y.Aju{2)Mju{l) 



(A4) 



where u{2) and m(1) are Dirac spinors, A's the invariant amplitudes and M's the gauge and 
Lorentz invariant matrices defined by 

Ml = -75^y (A5) 

M2 = 2-f5{e ■ pik ■ p2 - e ■ p2k ■ pi) (A6) 

M3 =75(A-Pi- Ve-Pi) (A7) 

M4 = 75(^-^2 -^e-P2) (A8) 



The F's are only defined in the CMS frame, but the A's are valid in an arbitrary frame. 
One can find the relations between the F's and the A's by first expanding Eq. ( |A^ ) in Pauli 
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space, then expressing the result in the CMS frame and putting it in the form of Eq. ( [A2D . 
After some algebra, we find 



_ /' E2 + 1712 

^ ~ 47r V 2W 
^ ~ An 



2W 



kq f E2 + 1712 



W + mi , 
-cAi + CA3 

W - mi , 
cAi + cAg - 



u — rrin 



2{W -mi 



-cAa 



u — mr, 



2{W -mi 



-cAa 



2W 



) [{W - mi)cA2 - cA^\ 



(A9) 
(AlO) 

(All) 



where c = [Ei + mi)(i?2 + m2)/(4mim2). 
Inverting these equations, we obtain 



cAi 
CA2 
cA^ 
cAa 



kV2W 
Air 



W — mi W + mi 

Fi H — , ^ -F2 



^E2 + m2 
F3 



kqV2W 
An 



y/E2 - m2 

F, 



kV2W 
An 



^E2 + m2 \JE2 - m2 
2Fi 2F2 



+ 



^E2 + m2 \JE2 - m2 
-F3 



u-m\ ( V^2 - ^2 ^ , V^2 + m2 „ ^ 

--^3 + 



W - mi 



W + mi 



kqV2W 



{W + mi)y/E2 + m2 {W - mi)y/E2^^m2 



(A12) 
(A13) 

(A14) 
(A15) 



Substituting the A's back into Eq. (|A4|), we obtain a transition operator which is valid in 
any reference frame and which is suitable for use in nuclear calculations. Eq. ( [A^ ) can be 
further expanded into Pauli space and cast into the form of a spin non-flip term plus a spin 
flip term 



Tfi = L + i cr ■ K.. 



(A16) 
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FIG. 1. Total cross section for the process (7,??) on the proton calculated with PS and PV 
Born terms. The full curve contains no Born terms, while the dashed lines are (from the top down) 
obtained with g'^j^j^/ At: =0.1, 0.5, 1.0, and 3.0 for PS-coupling, and g'^j^j^/ 477=1.0, 3.0, 6.0 and 10.0 
for PV-coupling, respectively. The experimental data are from Mainz | pT[ . 
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FIG. 2. Differential cross section for eta photoproduction at 766 MeV photon energy. The sohd 
and dashed hnes are calculations in PS coupling with coupling constant 0.4 and in PV coupling 
with 10, respectively. The data are from Mainz [pTj. 
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FIG. 3. Complete view of the differential cross section and the photon asymmetry of eta pho- 
toproduction on the proton for PS coupling with g^j^pf/Air = 0.4. 
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FIG. 4. The coordinate system for the quasifree reaction ^(7, r]N)B 
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FIG. 5. The energy dependence of the r] optical potential parameter b for two different models 



of the rjN t-matrix. See Sec. pID| for the discussion of DWl and DW2. 
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FIG. 6. Effects of distortions in the coincidence cross section and the photon asymmetry for 
^(TjW); P3/2 at £^=750 MeV, 9jf = 20° and 9p = 15°. The dashed hne is the plane wave 
calculation, the solid line is with ij distorted by the model DWl while the proton is kept as plane 
wave, and the dot-dashed line is with p distorted while rj is kept as plane wave. 
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FIG. 7. Effects of the r/-nucleus potential as a function of eta angle for plane wave, full DWl, 
and DWl with real part turned off. 
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FIG. 8. Effects of medium modifications to the S'ii(1535) resonance are shown by changing its 
bare mass by 3%. The reaction is ^'^C{'y,r]py^ Bg,s. at r^=100 MeV, 6*^ = 20° and Op = 20°. 
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FIG. 9. Same as in Fig. |, except for L>i3(1520). 
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FIG. 10. Eta energy dependence of the inclusive cross section (f a / dTj^dO^r^ on ^^C and "^^Ca 
calculated in the plane wave (PW) and the two different 77-nucleus potentials DWl and DW2. 
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FIG. 12. Comparison of our calculations with the Mainz data [12| for the inclusive cross section 
da/dT^ on 12C7 and ^°Ca. 
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FIG. 13. Comparison of our calculations with the Mainz data [12| for the inclusive cross section 
da/de^ on and ^^Ca. 
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FIG. 14. Comparison of our calculations with the Mainz data for the total cross section on C 
and 40Ca. 
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